Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 15 June 2009 (MN WI^ style file v2.2) 



Cluster strong lensing in the Millennium Simulation: The 
effect of galaxies and structures along the line-of-sight 



Ewald Puchwein^ and Stefan Hilbert^'^ 



o 
o 

in 



O 

u 

o 

(N 
> 
cn 
in 

(N 

^' 

O 

o\ 
o 



^ Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild-Strafie 1, 85741 Garching, Germany 
Argelander-Institut fiir Astronomie, Auf dem Hiigel 71, 53121 Bonn, Germany 



15 June 2009 



ABSTRACT 

We use ray-tracing through the Millennium simulation to study how secondary matter 
structures along the line-of-sight and the stellar mass in galaxies affect strong cluster 
lensing, in particular the cross-section for giant arcs. Furthermore, we investigate the 
distribution of the cluster Einstein radii and the radial distribution of giant arcs. We 
find that additional structures along the line-of-sight increase the strong-lensing optical 
depth by ~ 10 — 25%, while strong-lensing cross-sections of individual clusters are 
frequently boosted by as much as ^ 50%. The enhancement is mainly due to structures 
that are not correlated with the lens. Cluster galaxies increase the strong-lensing 
optical depth by up to a factor of 2, while interloping galaxies are not significant. 
We conclude that these effects need to be taken into account for predictions of the 
giant arc abundance, but they are not large enough to fully account for the reported 
discrepancy between predicted and observed abundances. Furthermore, we find that 
Einstein radii defined via the area enclosed by the critical curve are 10 — 30% larger 
than those defined via radial surface mass density profiles. The contributions of radial 
and tangential arcs to the radial distribution of arcs can be clearly distinguished. The 
radial distribution of tangential arcs is very broad and extends out to several Einstein 
radii. Thus, individual arcs are not well suited for constraining Einstein radii. 

Key words: gravitational lensing - dark matter - large-scale structure of the Universe 
- galaxies: clusters: general - cosmology: theory - methods: numerical 



1 INTRODUCTION 

In studies of strong gravitational lensing, the mass distri- 
bution deflecting the light is typically assumed to be con- 
centrated in a single region that is much smaller than the 
distance between source and observer. The extension of the 
matter distribution along the line-of-sight can be neglected 
in this case. The light deflection may then be described by 
projecting the lens' mass distribution onto a single plane and 
by assuming that light is deflected only when it passes this 
lens-plane and propagates along straight lines otherwise. In 
this so-called thin-lens approximation, the deflection angle 
is only sensitive to the projected surface mass distribution 
of the lens. 

This approximation has been widely used to study many 
aspects of gravitational lensing. In particular, it was adopted 
to investigate the efficiency of galaxy clusters to form giant 
arcs. However, studies employing it together with spheri- 



cally symmetric analytic cluster models (e.g. Wu and Mag 
1 19961 : iHattori et all 1 19971 : iMolikawa et~al] 1 19991 ) found sig- 
nificantly lower numbers of giant arcs than observed. The 
discre pancy is reduced when considering triaxia lly shaped 
halos l|Oguri et al ] |2003l : iMeneghetti et al.ll2007l ). More re- 



alistic lens models that also contain substructure can be 
obtained from cosmological simulations. They were used by 
several authors to i nvestigate the strong-lensing efficiency 
of clusters (see e.g. iBartelmann and Weiss! Il994l ). One of 
the findings was that even when using such detailed nu- 
merical models, the number of observed giant arcs is sig- 
nificant ljH^jge£_th^B_th*TOf^ in ^ ACDM cos- 
mology (|Bartelmann et ah I ll998l ). The discrepancy can be 
reduced by t aking cluster mergers llTorri et al.l 2004 ), bary- 
onic physics llPuchwein et al.ll2005l : IWambsganss et al]l2008l : 
iRozo et al.' '2008^) and a broad source redshift distribution 
(iWambsganss et al.u200'i : iLi et al.ll2005, ) into account. How- 
ever, the question whether or not the number of observed 
giant arcs is compatible with a ACDM cosmology is still not 
settled, especially for a low value of the m atter power spec- 
trum normalisation erg ^Fedeli et al." 2008) as suggested by 
recent cosmological constraints ()Komatsu et al...2009, ). 

Almost all of the studies mentioned above employed the 
thin-lens approximation and considered galaxy clusters as 
isolated gravitational lenses. In reality, however, clusters are 
embedded in a cosmic web of large-scale structures. In addi- 
tion, other collapsed objects like galaxies or galaxy groups 
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may be between a lensed background galaxy and the cluster 
or between the cluster and the observer and thus contribute 
to the lensing magnification and distortion. The effect of 
such additional structures along the line-of-sight has been 
ignored in most cluster- lensing studies. 

A notable exception, in which the impact of such 
secondary structures w as specifically investigated, is 
IWambsganss et ahl l|2005l ). Using a large cosmological simu- 
lation and multiple lens planes between source and observer, 
it was found that for high source redshift a significant frac- 
tion of the lines-of-sight with a convergence exceeding the 
critical value of unity for strong lensing/multiple image sys- 
tems, does so only due to a contribution from a second sub- 
dominant lens-plane. This suggests that additional struc- 
tures along the line-of-sight can, indeed, be important for 
strong lensing of high-redshift sources. From these results, 
it is, however, difficult to quantify the impact of such sec- 
ondary structures on t he abundance of giant arcs. 

On the other hand. lHennawi et al.l l|2007l ) also used mul- 
tiple lens planes in their analysis, but found that for sources 
at redshift Zs = 2, additional structures along the line-of- 
sight only mildly 7%) affect the number of multiple- 
image systems with image separations larger than 10 arcsec. 
But then, the largest effect was foun d for even higher source 
redshift in lWambsganss et al.l (|2005l ) and there may be a sig- 
nificant contribution from smaller-separation strong-lensing 
systems. 

There were also several studies investigating the im- 
pact of_clustcr galaxies on the lensing properties of clus- 
ters. [Meneghctti ot al. (2000) and Florcs ot al. (2000) found 
that the lensing efficiency of clusters is only mildly affected 
by its member galaxies. Brightest cluster galaxies, however, 
were not specifically accounted for in these studies, but were 
sh own to be able to raise lensing cross-sections by ~ 50% 
in lMeneghetti et al.l ((2003). Dalai ct al. (2004) showed that 
while the effect of galaxies on wide-separation giant arcs is 
small , it becomes s ignifica nt for smaller < 15" separation 
arcs. iHilbert et al.l (|2008lj) found that for point-like circu- 
lar sources, galaxies can boost the number of images with 
length-to-width ratios exceeding 10 by up to a factor of 2 
even when considering only images with separations > 5". 

In this work we investigate in detail how strong cluster 
lensing is affected by secondary matter structures along the 
line-of-sight and by both cluster and interloping galaxies. We 
carefully quantify their impact and explore its dependence 
on cluster mass, cluster redshift, and source redshift. We also 
study how the total lensing efficiency of a realistic cluster 
population is affected. Our strong-lensing simulations em- 
ploy ray-tracing along a backward light cone constructed 
from the Millennium simulation and use elliptical sources 
with sizes matching observations. The impact of galaxies is 
studied using a semi-analytic catalog that provides realistic 
galaxy properties and positions. 

We also compare different definitions of the Ein- 
stein radius, as it is a practical quantity for charac- 
terizing the strength of a gravitational lens. Recently, 
lOguri and BlandfordI (|2009l ) suggested to use the statistics 
of large Einstein radii as a cosmological probe, e.g. for mea- 
suring as. Here, we present the Einstein radius distribution 
of the Millennium simulation clusters which can be used 
to calibrate such methods. Finally, we investigate at what 
cluster-centric radii giant arcs are most likely to be found. 



The ray-tracing and strong-lensing simulation codes are 
introduced in Sect. [2] Our results are presented in Sect. [3] 
and summarized in Sect. |4] 



2 METHODS 
2.1 Ray-tracing 

Wo e mploy ray -tracing through the Millennium Simulation 
( Sp ringel et al.i 2005), a large Ai'-body simulation of cos- 
mic structure formation, to study the light defiection by 
matter structures in the observer's backward light cone. 
The Millennium Simulation assumes a fiat ACDM cosmol- 
ogy with a matter density flu = 0.25 (in units of the 
critical density), a cosmological constant with energy den- 
sity flA = 0.75, a Hubble constant of h — 0.73 (in units 
of 100 kms~^Mpc~^), and a scale-invariant initial density 
power spectrum with normalisati on gg = 0.9. The simu- 
lation used a TreePM algorithm (jSpringelllioOsI ) with 10^'' 
particles of mass nip — 8.6 x lO*/i~^M0 and a force soften- 
ing of 5h~^kpc comoving to simulate the structure formation 
in a cubic region of L = 500/i~^ Mpc comoving side length. 
Simulation snapshots were stored on disk at 64 output times 
between redshift z = 127 and z = 0. 

We use the mult iple-lens-plane algorithm described in 
iHilbert et all l|2008ah to calculate the gravitational defiec- 
tion of light by matter inhomogeneities between the source 
and the observer. The observer's backward light cone is di- 
vided into redshift slices. The matter content of each slice 
is then projected onto a lens plane transverse to the line- 
of-sight. Light rays are traced back from the observer to 
their source under the assumption that photons propagate 
unperturbed between these lens planes, but photons passing 
through a plane are deflected by an amount that is deter- 
mined by the projected matter distribution on the plane. 

The distribution of the dark matter in the observer's 
light cone is generated directly from the particle data of the 
Millennium Simulation. We use one redshift slice for each 
simulation snapshot. Each slice is filled with the simulation 
particles of the corresponding snapshot. The particles are 
projected onto a hierarchy of meshes on the correspond- 
ing lens plane. The matter distribution on the meshes is 
smoothed by an adaptive scheme that retains a resolution 
of about 5ft~^ kpc comoving in dense regions (e.g. in the cen- 
tres of massive dark-matter halos). The light deflection due 
to the dark matter is then calculated from the smoothed dis- 
tribution with Fast-Fourier- Transform methods, flnite differ- 
encing, and bi-linear interpolation. 

The distribution of the luminous matter in the light 
cone is i nferred from t he semi-analytic galaxy-formation 
model bv lOe Lucia and B laizot (2007) . The galaxies within 
each redshift slice are projected onto the same lens plane 
as the simulation particles. For each galaxy, the projected 
stell ar matter is assumed to follow a de-Vaucouleurs pro- 
file (|de Vaucouleurslll94^ 1 for the bulge component, and an 
exponential surface density profile for the disc component. 
The light defiections induced by stars in the galaxies are 
then calculated by analytic expressions^ 

^ Details about the tre atment of the l i ght def iection by the stellar 
matter can be found in IHilbert et al.l fcoOSbt) . 
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In order to study the lensing properties of galaxy groups 
and clusters, we select a sample of massive dark-matter ha- 
los at three lens redshifts = 0.28, 0.62, and 0.99 in the 
Millennium Simulation. For each halo, we set up light rays 
starting from the observer with directions on a regular grid 
of 1024 X 1024 rays in a 6' x 6' field-of-view. The observer 
position is chosen such that the halo centre is contained 
in the light cone spanned by the rays. The rays are traced 
back through the series of lens planes, and the ray positions 
on each lens plane, which also serve as source planes, are 
recorded. These ray positions are then used to create lensed 
galaxy images as described in Sect. 12.21 

In this work, we want to investigate the influence of the 
stellar mass in galaxies on the strong-lensing properties of 
groups and clusters. We thus perform the ray-tracing not 
only as described above, but also by ignoring the contribu- 
tion from the stellar matter in galaxies to the light deflection. 
We then compare the results we obtain by either including 
or ignoring light deflection by the luminous matter. 

Furthermore we want to quantify the contribution from 
additional matter structures along the line-of-sight. There- 
fore in addition to the complete multiple-lens-plane ray- 
tracing (hereafter full ray-tracmg), we also perform ray- 
tracing in the following ways: In one set of ray-tracing simu- 
lations, we only take into account the light deflection by the 
lens plane containing the selected dark-matter halo {single- 
plane ray-tracing). These ray-tracing simulations still in- 
clude the effects of correlated structures close to the lens- 
ing cluster but neglect contributions from matter structures 
at different redshifts. In addition we perform another set 
of simulations, in which only the matter in a cubic region 
of 4/i~^ Mpc comoving side length around the halo centre 
is taken into account in the ray-tracing {cluster-only ray- 
tracing). 



2.2 Strong-lensing simulations 

We investigate strong lensing in these simulated light cones 
by placing sources on six different source planes at redshifts 
zs = 1.0, 1.5, 2.1, 3.1, 4.2, and 5.7. We then find the images 
of each source using the ray positions from the ray-tracing 
described above and calculate the length L, width W and 
magnification for each image. This allows a statistical char- 
acterisation of the strong-lensing properties of the matter 
structures in each light cone. 

For finding the images of a large number of 
source galaxies , we f ollow the method introduced by 
iMiralda-Escude (Il993a ^ and adapt ed to non-analytic mod- 
els by Bartelmann and WeissI l| 19941 ) and iBartelmann et al.l 
lll995h. Th e algorithm is described in some detail in 



IPuch wein et al.' ('2005'). A previous version is also discussed 
inlMcnc ghotti et al . ( 200Q ). 

First, the critical curves and caustics for each lens are 
determined by checking where the determinant of the Jaco- 
bian of the lens mapping changes its sign. Then, the algo- 
rithm adaptively places source galaxies on the source plane 
such that regions near caustics are sampled with a higher 
source density than regions far from caustics. A statistical 
weight for each source accounts for the varying source den- 
sity. Source galaxies are modelled as ellipses with equivalent 
radii of 0.85, 0.65, 0.45, 0.3, 0.25, and 0.25 arcsec at redshifts 



Zs = 1.0, 1.5, 2.1, 3.1, 4.2, and 5.7 respectivelyQ Source el- 
lipticities are assigned by randomly drawing ratios of major 
to minor axis from the interval [1,2]. 

Each of the ray-traced fields-of-view is covered by a grid 
of 4096 X 4096 pixels in the image plane. The source-plane 
positions of these pixels are calculated by bi-linear interpo- 
lation between the source-plane positions of the 1024 x 1024 
rays that were traced through the simulation as described 
in Sect. 12.11 All image pixels of a source are then found 
by checking which of the grid pixels correspond to a source- 
plane position that is enclosed by the ellipse representing the 
considered source galaxy. The image pixels are then grouped 
into images. 

The properties of an image are obtained by finding the 
image pixel pi that falls closest to the source centre when 
mapped to the source plane, the pixel p2 that is the far- 
thest from pi, and the pixel pa that is the farthest from p2. 
We then fit a circle through these three points and use the 
arc length from p2 to pa as the length of the image. The 
image area is calculated directly from the number of image 
pixels, while its perimeter is obtained by walking along the 
boundary pixels and summing up their mutual distances. 
We apply the same resolution corre c tions to image length 
and perimeter as in iPuchwein et~all (|2005l ). although they 
are less irnportant here due to the better image plane grid 
resolution|j We then compute a simple geometric figure (el- 
lipse, circle, rectangle or ring) with equal area and length 
to determine the image width W, which is approximated 
by the minor axis of the ellipse, the diameter of the circle, 
the smaller side of the rectangle or the width of the ring, 
respectively. We choose the type of the figure by comparing 
its circumference to the perimeter of the image. 

Finally, for all light cones and each source redshift, we 
quantify the efficiency to produce arcs, which we define as 
images with a length-to-width ratio L/W > 7.5. More pre- 
cisely, we determine strong-lensing cross-sections o-x,/i/k>7.5 
by summing up the statistical weights of the sources that 
have images with L/W > 7.5 and calculating the corre- 
sponding area in the source plane. If there is more than one 
arc for a source, we multiply the statistical weight of this 
source by the number of arcs. 



3 RESULTS 

In order to quantify the impact of additional matter struc- 
tures along the line-of-sight on strong cluster lensing, we 
selected in total 100 cluster-sized halos, only by their mass, 
from the Millennium simulation. Table [T] summarizes their 
distribution in redshift Zh and halo mass MJqo Q For each 



^ These source radii agree well with measured half-light radii of 
galaxies at these redshifts (Ferguson et al. 2004). 
^ We also performed several strong lensing simulations with an 
increased 8192 X 8192 grid resolution to check that results are ro- 
bust even for the smallest source sizes. In all cases, strong lensing 
cross sections changed by less than 10% compared to the lower 
resolution results. Also note that any small residual resolution 
dependence will affect the results obtained by different kinds of 
ray-tracing in the same way. 

* Throughout this paper, we define cluster mass as the mass 
within a spherical region with a mean density 200 times the crit- 
ical density of the universe at the cluster redshift. 
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mass bin 
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at 2L=0.28 
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at zl=0.62 


10 


10 


10 


4 





at 2l=0.99 


10 


10 


7 









Table 1. Mass range and number of studied clusters within the 
mass bins at redshifts = 0.28, 0.62, and 0.99. We aim to use 
10 clusters per bin, but there are < 10 clusters in the Millennium 
simulation in the mass range of some bins. Also within each bin, 
clusters are roughly constantly spaced in logarithmic cluster mass. 

cluster, we then performed full, single-plane, and cluster- 
only ray-tracing simulations, as described in Sect. 12.11 and 
calculated the corresponding strong-lensing cross-sections as 
detailed in Sect. \2l2\ Comparing them to each other allows 
us to clearly pin down the effect of additional line-of-sight 
structures on cluster lensing properties and giant arc statis- 
tics. 

For all selected clusters, we also performed full and 
cluster-only ray-tracing simulations that ignored the light 
deflection by the stellar matter in both cluster and interlop- 
ing galaxies. This enables us to also quantify the impact of 
the stellar component of galaxies on strong cluster lensing. 

3.1 The impact of additional structures along the 
line-of-sight 

In Fig. [l] we show the average strong-lensing cross-sections 
f"L/H'>7.5 for all cluster mass bins given in Table [l] and all 
considered source and lens redshifts. For each mass bin and 
lens redshift, a set of three curves illustrates the results ob- 
tained when using the full, single-plane, and cluster-only 
ray-tracing. In all cases, the stellar matter in galaxies was ac- 
counted for in the ray-tracing. In the left panels, all suitable 
arcs in the whole 6' x 6' ray-tracing fields where considered 
for calculating cri/vi/>7.5. 

Probably the most obvious result is that the average 
cross-section is a strongly increasing function of source red- 
shift for all mass bins, lens redshifts, and types of ray- 
tracing. Clusters are, thus, more efficient in lensing source 
galaxies at higher redshi ft. This finding is in good agreement 
with previous studies llWambsganss et al.l 120041 : Li et al.l 
l2005l : lllilbert et al.ll2007l l. 

A more careful inspection of the results reveals that for 
all considered lens redshifts, the cross-sections of the most 
massive clusters (uppermost set of curves in each panel) are 
only very mildly affected by additional matter structures 
along the line-of-sight. The results obtained using full ray- 
tracing are very similar to those found by single-plane and 
cluster-only ray-tracing for such massive clusters. 

The effects of additional matter structures along the 
line-of-sight are more important for less massive clusters 
and for sources at high redshift. There, the cross-sections 
obtained when using the complete multiple-lens plane ray- 
tracing are on average ~ 20 — 40% larger than those found 
by considering only light defiection by the cluster. On the 
other hand, the difference between results from single-plane 
and cluster-only ray-tracing are still very small. Thus the 
increase in lensing efficiency found when accounting for all 



the matter in the light cone is not due to correlated struc- 
tures near the cluster, but due to uncorrelated structures 
along the line-of-sight. 

So far, all suitable arcs, i.e. with a length-to- width ratio 
exceeding 7.5, found in the ray-traced field were counted for 
the computation of a cluster's strong-lensing cross-section. 
However, an observer might not assign all these arcs to the 
selected cluster, but also to other objects that happen to be 
in the ray-traced field. Neglecting this is unlikely to signif- 
icantly bias results for the most massive clusters. However, 
the smallest clusters in our sample are in some cases not 
even the most prominent lens in the field-of-view. Thus by 
assigning all found arcs to such a cluster, one would over- 
predict its lensing efficiency, especially when using the full 
ray-tracing. 

To prevent such misassignments from biasing our re- 
sults, we repeated the computation of o"_l/w>7.5 for all clus- 
ter, but this time counting only arcs within 5 Einstein radii 
from the selected cluster's centre. Note that here, we define 
the Einstein radius Oe hy 

Oe = \f^, (1) 

where A is the solid angle enclosed by the cluster's outer 
critical curve and calculated using cluster- only ray-tracing. 
In this way, we largely avoid assignments to a wrong object. 
In rare cases of ongoing mergers or very elliptic clusters with 
substructure, we might, however, miss some arcs that are 
very distant from the cluster centre but should probably 
be assigned to it. We employ this cut in image separation 
troughout the remainder of this paper. 

Results obtained in this way are shown in the right pan- 
els of Fig. [T] As expected, the cross-sections of the most 
massive clusters are almost unchanged. For the smallest clus- 
ters on the other hand, the lensing cross-section obtained by 
full ray-tracing is reduced and becomes more similar to the 
cluster-only results. Thus, when counting only arcs within 
5 Einstein radii, most small clusters are not significantly af- 
fected by additional structures along the line-of-sight. The 
effect is biggest for intermediate mass clusters, where cross- 
sections increase on average by roughly ~ 15 — 30%. 

For individual clusters, however, effects can be signifi- 
cantly larger. We illustrate this in Fig. [21 which shows the 
ratio of the cross-section obtained by the full ray-tracing to 
the values found by cluster-only ray-tracing. For the sake of 
clarity, only the results for the 2l=0.28 clusters and three 
different source redshifts are shown. Again, only arcs within 
5 Einstein radii were counted. Boosts of the cross-sections 
by additional matter structures of up to 50% happen fre- 
quently. There are also a couple of clusters with significantly 
larger ratios. For most of them, however, absolute values of 
the cross-sections are quite small, so that a small absolute 
increase results in a large relative boost. 

A particularly impressive example is the second most 
massive cluster shown in the figure. In spite of its enormous 
mass, secondary structures along the line-of-sight boost its 
cross-section for lensing high-redshift sources by roughly 
60%. A large part of this increase is due to three group-sized 
halos, two at redshift 1.5 and one at redshift 1.63, which are 
close to the cluster in projection. In addition, there is a sig- 
nificant effect due to a 2 x IO^^Mq/H halo that is also at 
redshift 1.63. Normally, one would not expect such a small 
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Figure 1. Average strong-lensing cross-sections crL/w>7.5 of clusters in mass bins (a)-(e) (see Table[T}. The cross-sections were obtained 
using full light cones (red solid lines), only one lens plane containing the selected cluster {green long-dashed), and light deflection only 
by the selected cluster {blue short- dashed). Results for clusters at redshift 2l=0.28 {bottom panels), zj^ = 0.62 {middle panels), and 
2l = 0.99 {top panels) are shown. In the left panels, arcs in the whole 6' X 6' field were considered, while in the right panels, only arcs 
within 5 Einstein radii were included. 



object to have a noticeable impact on the lensing properties 
of a cluster whose mass is almost three orders of magnitude 
larger. This halo is, however, itself strongly lensed by the 
foreground galaxy cluster. The peak in the lensing conver- 
gence corresponding to it is shaped like a long thin tangential 
arc. Due to the lensing magnification of this halo, its light 
deflection affects a much larger number of lines-of-sight. In 
this way, this small object can si gniflcantly alt e r the c riti- 
cal curve of the cluster. Also see iHilbert et al. I (|2008ah for 
a discussion of how multiple-lens-plane ray-tracing can be 
used to study lens-lens coupling effects in weak lensing. Fi- 
nally, there is also some contribution from several halos with 
masses of a few 10^^Mq//i and redshifts between zl =1.7 
and 3.3. Remarkably, the importance of additional line-of- 
sight structures is already strongly reduced for sources at 



a somewhat lower redshift of Zs=3.1. This is because the 
efficiencies of lenses at zl > 1-5 significantly increase with 
source redshift between 2s=3.1 and 5.7, while there is only 
a marginal gain for a lens at the much lower cluster redshift 
of ZL =0.28. 

There are also several clusters which have a smaller 
cross-section when including the effects of additional struc- 
tures along the line-of-sight. However, the decrease is less 
than 10% in most cases. Note that such clusters are not 
necessarily projected onto underdense regions, as we do not 
find a clear correlation between change in cross-section and 
change in average convergence within 5 Einstein radii0 In 



The region inside which the average convergence was computed 
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Figure 2. Ratio of the strong-lensing cross-section o"i,/w>7,5 
found when using the full multiple-lens-plane ray-tracing to the 
value obtained when considering only light deflection by the se- 
lected cluster. Results are shown as a function of cluster mass 
for all zl=0.28 clusters in our sample. Three different source 
redshifts, zg = 1.0, 3.1, and 5.7 were used. Only arcs within 5 
Einstein radii were included in the analysis. 



most cases where we find significantly smaller cross-sections 
when including all structures in the light cone, it is the 'ex- 
ternal' shear caused by the secondary structures along the 
line-of-sight that counteracts the cluster shear on potential 
arcs. In many cases, this lowers their length-to-width ratios 
below the adopted threshold value of 7.5. 

Figure |3] shows the critical curves and caustics of two 
clusters and illustrates how they change when the effects 
of additional structures along the line-of-sight are included. 
The first cluster, shown in the figure's upper panels, is an ex- 
ample of a strongly affected object. Its lensing cross-section 
is boosted by a factor of two, and the shape of its criti- 
cal curve changes significantly. The latter is stretched and 
merges with the critical curve of a subhalo. Also shown are 
the centres of all giant arcs formed in the cluster's full and 
cluster-only strong-lensing simulations as well as the sources 
that give rise to them. They illustrate how the shape of the 
cluster's strong-lensing cross-section changes and where ad- 
ditional arcs appear. The second cluster is a more typical 
example. Its cross-section increases by about ~ 30%, mostly 
due to an additional structure below and to the right of 
the cluster. The shape of its critical curve and caustic are 
only mildly affected. One can however clearly see that the 
caustics are shifted with respect to each other due to light 
defiection by structures between the cluster and the source 
plane. 

As mentioned above and illustrated in Fig. [1] the im- 
pact of additional structures along the line-of-sight on clus- 
ter strong-lensing cross-sections depends on cluster mass. 
Thus, one needs to take the steep cluster mass function into 
account when investigating how the total lensing efficiency 
of a realistic cluster population is affected. However, instead 
of employing an analytic mass function for that purpose, we 



was kept fixed, as we used the Einstein radius obtained from 
cluster-only ray-tracing in both cases. 



directly use the full Millennium simulation halo catalogue. 
We start from the cross-sections we calculated for our se- 
lected cluster sample and then assign a cross-section to all 
other Millennium simulation clusters, here defined as halos 
with M200 > W^'^h~^ M0, by interpolating linearly in clus- 
ter mass. Summing over all clusters finally yields an estimate 
of the total strong-lensing cross-section of the Millennium 
simulation cluster population at zl = 0.28, 0.62, and 0.99. 
We calculated these quantities based both on full and on 
cluster-only ray-tracing. Figure |3] illustrates how additional 
structures along the line-of-sight affect the lensing efficiency 
of the cluster population as a whole and thus the contribu- 
tion of clusters at these redshifts to the strong-lensing op- 
tical depth. Clearly, the lensing efficiency for high-redshift 
sources increases when taking such additional structures into 
account. This effect is larger for intermediate and high red- 
shift cluster populations, where we find boosts up to 30%. 
For the = 0.28 clusters and high redshift sources, the 
lensing efficiency increases by about 10%. 

We also show 68% confidence intervals of the total cross- 
section ratios. For each cluster redshift, they were calculated 
from 10^ bootstrap cluster samples. Note that we bootstrap 
both the clusters which we use for the cross-section inter- 
polation as well as the full cluster population in the simula- 
tion. The former accounts for uncertainties due to the lim- 
ited number of objects for which we have performed strong- 
lensing simulations and for noise in the cross-section calcu- 
lation, while the latter accounts for cosmic variance in the 
whole simulationlfl Most uncertain is the contribution from 
the large number of low mass clusters to the strong-lensing 
optical depth. This is because their cross-sections are small, 
especially for low source redshift, and thus the relative error 
in their determination is larger. 

Of course, the strong-lensing cross-section of individ- 
ual clusters does not only depend on cluster mass, but due 
to halo triaxiality, also on orientation (iDalal et al.li2003). 
as w ell as on halo concent ration (see e.g. 



Hennawi et al.l 



l2007t ). and dynamical state (|Torri et al.ll2004 ). It is, there 
fore, not suprising that the distribution of cross-sections is 
quiet broad even at constant cluster mass. Obviously, the 
impact of structures along the line-of-sight also depends on 
the line-of-sight 's orientation and can change significantly 
depending on it. Thus, while the curves given for individual 
mass bins in Fig. [T] (and Fig. [S] in the next section) are very 
useful for identifying trends with cluster mass, they may, 
due to the limited number (^ 10) of objects per bin, not 
always be fully representative for all cluster lenses in that 
mass range. The results we give for the whole cluster pop- 
ulation, on the other hand, should be much more robust as 
a larger number of clusters is used for their derivation. Also 
note that the bootstrap error estimates include the effects 
mentioned above. 



^ Before we applied bootstrapping, we have performed a Monte- 
Carlo analysis that employed random cross-sections drawn based 
on a rough estimate of the cross-section cluster mass relation and 
its scatter to confirm that bootstrapping gives reliable error esti- 
mates even when using such an interpolation scheme to calculate 
total cross-section ratios. 
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Figure 3. Critical curves {left panels) and caustics (right panels) of a 10^^/i~^Mq cluster at 2l = 0.62 (upper panels) and a 7 X 
lO^*/!"^ Mq cluster at zj^ = 0.28 (lower panels). Curves obtained by considering light-deflection only by the cluster (blue) as well as by 
all matter structures along the line-of-sight (red) are shown. A source redshift of zg = 3.1 was assumed. The cluster shown in the upper 
panels is a prominent example for an object whose lensing efficiency is strongly boosted by additional structures along the line-of-sight. 
Its cross-section increases by almost a factor of 2. The lower panels show a more typical example, where the cross-section increases by 
~ 30%. In both cases only giant arcs within 5 Einstein radii (green circles) were counted. In the left panels, the centres (defined by 
pixel pi as described in Sect. 12.21 1 of all arcs with L/W > 7.5 found in the strong-lensing simulation are shown. In the right panels, the 
positions of the sources that give rise to these arcs are indicated. Note that for the sake of clarity, we choose a dilferent scale in the 
bottom, right panel. 



3.2 The impact of the stellar mass in galaxies 

In Fig. [H we compare the strong-lensing cross-sections cal- 
culated accounting for the light deflection by the stellar com- 
ponent of galaxies to those found ignoring all stellar mat- 
ter. Average cross-sections for all the zl ~ 0.28, 0.62 and 
0.99 cluster mass bins listed in Table [1] are shown. Results 
obtained by both full and cluster-only ray-tracing are indi- 
cated. When accounting for stellar matter, the former in- 



clude the effects of all galaxies between source plane and 
observer while the latter include only the effects of cluster 
galaxies. One can clearly see that the stellar matter signifi- 
cantly increases the strong-lensing cross-sections of clusters, 
especially for low-mass systems and for low source redshift. 
In other words, the effect is largest for the most inefficient 
lenses, while very massive clusters are only mildly affected. 
It is also worth noting that the impact of stellar matter is 
very similar when using full and cluster-only ray-tracing. 
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Figure 4. Impact of additional structures along the line-of-sight 
on the combined strong-lensing cross-section of all Millennium 
simulation clusters. Ratios of total cross-sections estimated from 
full and cluster-only ray-tracing are shown as a function of source 
redshift for the cluster populations at zl = 0.28, 0.62, and 0.99 
(thick lines). Only arcs within 5 Einstein radii were included in 
the analysis. Also shown are 68% confidence intervals of the ratios 
(thin lines) obtained by bootstrapping techniques. 



This means that most of the increase is due to the steUar 
component of cluster galaxies, while interloping galaxies be- 
tween source and cluster or between cluster and observer 
play only a minor role. 

Also in this context, the impact on cluster lensing 
cross-sections depends on cluster mass. We, thus, use the 
same method as in Sect. 13.11 to estimate how the lens- 
ing efficiency of a realistic cluster population is affected. 
Figure [6] illustrates the increase of the combined strong- 
lensing cross-section of all Millennium simulation clusters 
with M|Jo > 10^" /i"^ Mq. Relative changes due to the stel- 
lar components of galaxies are shown as a function of source 
redshift for all clusters at zl ~ 0.28, 0.62 and 0.99. We also 
plotted the 68% confidence intervals of the cross-section ra- 
tios obtained by a bootstrap error analysis similar to that 
described in Sect. 13. Il Full ray-tracing was used to derive all 
the data shown in the figure. Stellar matter boosts the lens- 
ing efficiency of the ~ 0.28 cluster population by almost 
a factor of two for low redshift sources. For high-redshift 
sources the strong-lensing efficiency still increases by about 
50%. The clusters at higher redshift show the same trend in 
source redshift with somewhat larger boosts for low redshift 
sources. However, the magnitude of the effect is only about 
half of that found for the zl = 0.28 clusters. 

In order to understand this difference better, we cal- 
culated the Einstein radius of the average strong-lensing 
cluster for all lens redshifts. More, precisely we assigned an 
Einstein radius to all Millennium simulation clusters using 
the same interpolation method as for the cross-sections. We 
then derived the cross-section-weighted Einstein radii for the 
three cluster populations and all source redshifts. We found 
that for sources at redshift 1.5 the Einstein radius of the av- 
erage lensing cluster is very similar for clusters at zl = 0.28 
and 0.62, namely about ~ 8 arcsec, and larger in the latter 
case for higher redshift sources. For the zl ~ 0.99 clusters 
it is compareable or larger than for the zl ~ 0.28 clusters 
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Figure 5. Average strong-lensing cross-sections o'L/wyT.S °f 
clusters in mass bins (a)-(e) at cluster redshift zl=0.28 (bottom 
panel), zl=0.62 (middle panel), and 2l=0.99 (top panel). The 
cross-sections were obtained using full ray-tracing and accounting 
for stellar matter (red solid lines), cluster-only ray-tracing and 
accounting for stellar matter (dark-blue, dashed), full ray-tracing 
and ignoring stellar matter (magenta, dotted), cluster-only ray- 
tracing and ignoring stellar matter (light-blue, dot-dashed). Only 
arcs within 5 Einstein radii were counted. The red solid and dark- 
blue dashed curves show the same data as in the middle right and 
lower right panels of Fig. [T] 



above a source redshift of 2.1. Thus, in spite of the much 
larger distance between lens and observer the Einstein angle 
of strong-lensing selected clusters does not decrease at high 
lens redshift. It corresponds, however, to a larger physical 
distance at the cluster position. The stellar mass within the 
critical curve is typically dominated by the central galaxy, 
which is more concentrated than the dark matter halo. Thus, 
when the physical size of the critical curve increases lensing 
becomes more dark matter dominated and less affected by 
stellar matter. 

Our results for the increase of the lensing cross sec- 
tions due to cluster galaxies are in good agreement with 
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Figure 6. Impact of stellar matter on the combined strong- 
lensing cross-section of all Millennium simulation clusters. Ratios 
of the total cross-sections estimated accounting for and ignoring 
stellar matter are shown as a function of source redshift for the 
cluster populations at = 0.28, 0.62, and 0.99 (thick lines). Also 
shown are 68% confidence intervals of the ratios {thin lines) ob- 
tained by bootstrapping techniques. The analysis employed full 
ray-tracing. Only arcs within 5 Einstein radii were included. 



2x10' 



the results of iHilbert et af] ()2008bl l for point-like sources. 
At first glance, these results may seem to be at odds with 
the findings of Mene ghetti et all (|2000l ). However, we find 
the largest effect in small clusters which were not consid- 
ered there. In addition, we also accounted for brightest clus- 
ter galaxies which were also not included in that study, but 
shown t o be able t o rais e lensing cross-sections by ^ 50% in 
[Meneghetti et all (|2003l ). 



3.3 Optical depths for giant arcs 

We can use our cluster sample to estimate the optical 
depth Ti/n/>7.5 for arcs with L/W > 7.5 as a function 
of source redshift zs- For the considered cluster lens red- 
shifts zl = 0.28, 0.62, and 0.99, we calculate the cross- 
sections o-i^/w>7.5,i of every halo in the simulation with mass 
M200 ^ IO^^/i'^Mq by linearly interpolating the cross- 
section as a function of logarithmic cluster mass between 
the values of the ray-traced clusters. 

We then compute an estimate of the differential opti- 
cal depth dTL/yY>7.5/dxL at the comoving line-of-sight dis- 
tances XL = x(2l) of these three redshifts 2l by 



'L/W>7.h 



2 X^i 0"i/W>7.5,i 

Xl- 



14o 



(2) 



where the o-L/w>7.5,i ^-r^ assumed to be in steradians 
and the summation goes over all clusters with M|oo ^ 
10^*/i~^ Mq found in the simulation output at redshift zl. 



1x10 




Figure 7. Differential optical depths 9t2^/w>7.5/9^l for arcs 
with length-to-width ratio L/W > 7.5 as a function of lens reshift 
2l for various source redshifts zg. The values obtained from full 
ray-tracing including the stellar matter and considering only arcs 
within 5 Einstein radii of the lens cluster are shown as symbols. 
The lines indicate the differential optical depths obtained from 
linear interpolation in the comoving line-of-sight distance xl- 



Q Vbox = {500h~^ Mpc)^ is the comoving volume of the sim- 
ulation box. 

Furthermore, we assume that 9tj^/w>7.5/9xl vanishes 
for zl = 0, zl ^ Zs, and zl > 2.6 (there are no clusters with 



> for ZL > 2.6). The contribution from 



other redshifts is then calculated by linear interpolation in 
XL- The resulting differential optical depths are illustrated 
for all considered source redshifts in Fig. [7] We show them as 
a function of the lens redshift zl to make them more easily 
comparable to other works. 

The optical depth Ti^^wyr.s is obtained by integrating 
drL/wyv.s/dxL between Xl = and xl = xi^s)- Simi- 
larly, we calculate the optical depth Ti^/wyio for Btics with 
L/W > 10. The optical depths resulting from full ray- 
tracing including the stellar matter in galaxies is compared 
to the estimates based on cluster-only ray-tracing with and 
without stellar matter in Fig. [S] 

The qualitative behaviour of t"£/m/>7.5 B.nd t^/vi'>io is 
very similar. The additional matter along the line-of-sight 
increases the optical depths by 10% for zs = 2.1 and by 
25% for zs = 5.7. Hence, projection effects are more im- 
portant for higher source redshifts. The stars in galaxies 
enhance the optical depths by a factor two for zs = 1 and 
cause an increase by 25% for sources at zs = 5.7. Hence, 
the stellar matter is more important for lower source red- 
shifts. Together, the stars and the additional matter along 
the line-of-sight increase the optical depth by a factor 1.5 to 
2. 

For comparison with iFedeli et all l|2008h . we compute 



Note that the contribution of low-mass clusters to the opti- 
cal depth is more important for high redshift sources. Our data 
suggest that for halos with M^qq < 10^*/i~^ Mq, it becomes non- 
negligible for zg ^2.1. This contribution is not included in this 
analysis. On the other hand, it is likely to strongly depend on 
the exact selection criteria for arcs. For example requiring some 
minimum value for the image separation would largely suppress 
it. 
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Figure 8. Optical deptlis Ti/vi/>7.5 ((a), top panel) and 
Ti^/lYyld {(h), bottom panel), for arcs with length-to-width ra- 
tio L/W > 7.5 and L/W > 10 as a function of source redshift 
zs ■ Shown are the values for full ray-tracing including the stellar 
matter {red solid lines), for cluster-only ray-tracing with stars 
{blue dotted lines), and for cluster-only ray-tracing without stars 
{green dashed lines). Only arcs within 5 Ein stein rad i i were in- 
cluded. Also shown is t^/vi/>io computed by I Li et al.l ||2005| ) for 
sources with 1" effective diameter {purple dash-dotted line). 

average optical depths 

t"l/W>7.5 ^ J 'izsPs{zs)TL/W>7.5(zs) (3) 

using the source redshift distribution (ISmail et al.lll995l ) 

ps{zs) = ^zs exp {~4^^) ■ (4) 

The resulting values Tj^/vi^>7.5 = 2.5 x 10^'' for full ray- 
tracing and fL/w>7.5 = 1-7 x 10~® for cluster-only ray- 
tracing withou t stars are 4 to 6 t imes smaller than the op- 
tical depths of iFedeli et al.l l|2008l ). Part of the discrepancy 
may be due to the fact that we assumed larger source radii 
at low redshifts. In addition, there may be some difference 
in the lensing properties between the simulated cluster pop- 
ulation used in this work and the analytic cluster popula- 
tion model employed there. In particular, cluster mergers 
are accounted for in different ways. In our simulations, clus- 
ter mergers naturally occu r. However, we mig ht be limited 
by small- number statistics. iFedeli et al.l (|2008| ). on the other 
hand, explicitly include cluster mergers, assuming that they 
proceed at the gravitational free-fall timescale and can be 
represented by two NFW halos approaching at constant ve- 



locity. Without the merger contribution their optical depths 
would be ~ 3 ti mes smaller and agree m uch better with our 
results. However. iHennawi et al.l (200'f) question that merg- 
ers can significantly boost optical depths for giant arcs. Con- 
sidering all these points, it seems likely that at least part of 
the discrepancy between our results and those of lFedeli et al.l 
(2008) is related to the way cluster mergers are accounted 
for, i. e. the importance of merger boosti ng of strong l ensing 
cross-sections may be overestimated in IFedeli et al.l (|2008h 
and/or our selected cluster sample may still be too small 
to fully account for strong boosts during rare major cluster 
mergers. 

In Fig.[SIb), we compare ou r optical depths for arcs with 
L/W > 10 to those obtained bv lLi et al.l (|2005h . The optical 
depths are very similar for sources at 2s = 1- However, our 
optical depths are 3 to 4 times larger for zs 55 2.1. The fact 
that we used smaller source radii at h igher redshifts accounts 
only for an increase of about 50% l|Li et ahllioosh . A rea- 
son for the remaining discrepancy could be the significantly 
better spatial resolution of the Millennium Simul ation com- 
pared to the simulations used by iLi et al.l (|2005t ) , which is 
likely to boost the lensing efficiencies of poor clusters, whose 
contribution to the optical depth is most important at high 
source redshift. Another reason may be the differences in 
the assumed cosmology. 

Admittedly, there is some uncertainty in our optical 
depth estimates due to the rather rough sampling of the 
cluster population, especially in cluster redshift. This may 
contribute to some of the discrepancies menti oned above. 
However, as our results deviate from those of IFedeli et al.l 
( 2008*) and lLi et al.l (|2005l ) in opposite directions, it is cer- 
tainly not the only cause. Also the relative changes in the op- 
tical depths due to including secondary structures along the 
line-of-sight and the stellar components of galaxies should be 
robust as sparse sampling will affect all results in a similar 
way. 

3.4 Einstein radii 

It is straightforward to define an Einstein radius for a spher- 
ically symmetric lens, but there are several ways to gener- 
alise the concept of an Einstein radius to realistic, i.e. non- 
spherical, cluster lenses. In Eq. ((TJ, we defined the Einstein 
radius 9e of a cluster via the solid angle enclosed by its 
outer critical curve. Another common definition is the ra- 
dius 6'^ of a circular region around the cluster centre having 
a mean projected surface mass density equal to the critical 
density for strong lensing. In the following, we compare the 
two definitions. 

We computed 9'^, which we conveniently obtain from 
each cluster's projected surface mass density profile, for the 
750 most massive clusters in the Millennium simulation. Fur- 
thermore we used single-plane ray tracing to calculate 9e 
for the 100 most massive Millennium halos and all additional 
lower mass clusters available in the cluster sample that we 
used in the previous sections. Comparing the results for all 
clusters for which both values were derived, we find that 

* We are most interested in large Einstein radii produced by 
massive clusters, which are not strongly affected by additional 
structures along the line-of-sight. 
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Figure 9. Distribution of Einstein radii of the 750 most mas- 
sive Millennium clusters at ^l=0-28 and sources at zg = 2.1. 
Distributions are given for the Einstein radius d'^, derived from 
surface density profiles (solid line), and 9e) derived from the crit- 
ical curves {dashed line). Also shown is the ^B-distribution of the 
clusters the SE-interpolation was based on {dotted line). 



the Einstein radii defined via the critical curve are on aver- 
age 10 — 30% larger (depending on cluster mass) than the 
radii obtained from the surface density profiles. For example, 
the mean and la-scatter of the ratios Se/^b for sources at 
23=2.1 and clusters at redshift 2;l=0.28 in mass bins (a)-(e) 
are 1.12±0.12, 1.15±0.18, 1.18±0.18, 1.32±0.11, and 1.28, 
respectively. The scatter is mostly uniform with only very 
few large outliers, which are related to merger processes. 

The difference between 6'^ and 0e can also be seen in 
Fig. |9l which shows the distribution of Einstein radii 6'^ and 
0E of the 750 most massive clusters in the Millennium sim- 
ulation at 2L=0.28, i.e. for all clusters at that redshift with 
■^200 > 1-47 X 10^*/i~^ Mq. The eE-distribution is based on 
120 clusters for which 6^ was calculated employing single- 
plane ray tracing. All other clusters were assigned a ^B-value 
by linear interpolation in cluster mass. Obviously, the 6^- 
distribution is shifted towards larger values compared to the 
^E-distribution. The Einstein radius distribution functions 
shown here should also be valuable for calibrating semi- 
analytical model predict ions of the Einstein radius distri- 
bution such as those in lOguri and BlandfordI l|2009D . who 
suggested they might be used as cosmological probes. 

3.5 Radial distribution of giant arcs 

Here, we study at what cluster-centric radii long thin arcs 
are most likely to be found. For a perfectly spherical sym- 
metric halo, the longest arcs are expected to appear close 
to the Einstein radius. However, for a realistic cluster, this 
is, in general, not the case. Even prominent giant arcs that 
consist of multiple merged images of a background galaxy 
and thus cross the tangential critical curve may be signifi- 
cantly outside the Einstein radius due to the critical curve's 
ellipticity. Arcs that do not consist of multiple images can be 
found at even larger radii. This can be seen in the left panels 

^ No scatter is given for mass bin (e) as it contains only one 
object. 



< 0.04 




Figure 10. Radial distribution of arcs with L/W > 7.5 for 
sources at redshift 2g = 2.1 and clusters in mass bins (a)-(e) at 
redshift zl = 0.28. The angular distance is given in units of the 
Einstein-radius 9'^, which is calculated for each cluster from a 
projected surface mass density profile. The vertical lines indicate 
the Einstein-radii derived from the areas enclosed by the clusters' 
critical curves. Full ray-tracing was used. 



of Fig. 131 which show that the centres of arcs with length to 
width ratios exceeding 7.5 are typically located outside the 
most distant parts of the tangential critical curve. 

We investigate the radial distribution of arcs in more 
detail by calculating cross-sections for arcs with L/W > 7.5 
as a function of radial distance in specific radial bins. In 
Fig. 1101 the results are shown for the five cluster mass bins 
listed in Table [1] In total 39 clusters at redshift = 0.28 
and full ray-tracing were used in this analysis. The cluster- 
centric angular distance is given in units of the Einstein 
radius 9'^. All curves were normalized by the total strong- 
lensing cross-section for arcs at all radii. Also indicated are 
the mean Einstein radii obtained from the critical curves of 
the clusters in each binE3 

One can clearly distinguish the contribution of radial 
and tangential arcs to the radial distribution of arcs. The 
former are typically found a.t 6 < 0.6 ^e- In contrast, the 
radial distribution of tangential arcs is very broad, and ex- 
tends from ~ 0.6 6'e out to ~ 4&b- Using a larger threshold 
value for the length-to-width ratio makes the distribution 
only slightly narrower. For L/W > 10, it still extends to 
values > 3 ^E • This shows that individual tangential arcs 
are not well suited for constraining Einstein radii, instead 
more detailed cluster mass models are needed. It is, thus, 
important to keep the broad distribution of tangential arcs 
in cluster-centric radius in mind when interpreting strong- 
lensing observations of galaxy clusters. 



4 SUMMARY AND CONCLUSIONS 

We have used ray-tracing through the Millennium simula- 
tion to study how secondary matter structures along the 



The indicated Se/^e values differ slightly from the numerical 
values quoted previously in the text as full ray-tracing and a 
smaller cluster sample were used here. 
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line-of-sight affect strong cluster lensing, in particular the 
cross-section for giant arcs. We also investigated the impact 
of the stellar components of both cluster and interloping 
galaxies on cluster lensing efficiencies. Furthermore, we ex- 
plored the distribution of the cluster Einstein radii and the 
radial distribution of giant arcs. 

We performed ray-tracing along backward light cones 
using a multiple-lens-plane algorithm that takes light de- 
flection by all matter structures between source and ob- 
server into account. We also performed single- lens-plane ray- 
tracing simulations that either considered only light deflec- 
tion by the cluster itself or by all the matter contained in 
the same lens plane as the cluster. Thus in the last case, the 
effects of correlated structures were included while indepen- 
dent matter structures at different redshifts were not. 

Comparing the results obtained by the different kinds 
of ray-tracing allows us to clearly pin down the impact of 
additional matter structures along the line-of-sight: 

• The strong-lensing efficiency of clusters increases when 
including the effects of additional structures along the line- 
of-sight. 

• The enhancement is larger for higher source redshifts. 

• The enhancement is mainly due to structures along the 
line-of-sight that are not correlated with the lens. 

• For individual clusters, boosts of the strong-lensing 
cross-section of up to 50% occur frequently for sources at 
redshift 3 and higher. 

• The strong cluster lensing optical depth increases by 
10% for zs ~ 2.1 and by 25% for 23 — 5.7 when including 
the effects of additional structures along the line-of-sight. 

Using the Millennium simulation's semi-analytic galaxy 
catalogue, we also investigated how the stellar components 
of galaxies affect the lensing properties of clusters. Compar- 
ing the results obtained when neglecting and including light 
deflection by the stellar components of cluster and interlop- 
ing galaxies, we find: 

• The strong-lensing properties of massive clusters are 
only mildly affected by stellar matter. 

• The lensing cross-sections of less massive clusters can 
be significantly boosted by cluster galaxies. 

• The enhancement due to the stellar matter is largest 
for low cluster and source redshift. 

• Interloping galaxies do not affect the cross-section sig- 
niflcantly. 

• Including stellar matter boosts the strong-lensing opti- 
cal depth by a factor of 2 for zs — 1 and by 25% for 23 = 5.7. 

Together, the stars and the additional matter along the 
line-of-sight increase the strong-lensing optical depth by a 
factor of 1.5 to 2, depending on source redshift. 

Furthermore, we computed the distribution of Einstein 
radii using two different deflnitions of the Einstein radius. 
We also determined the radial distribution of arcs compared 
to the Einstein radius of the lensing cluster. We obtain the 
following results: 

• The Einstein radii deflned via surface mass density pro- 
flies and via the area enclosed by the critical curve differ by 
10% to 30%. 

• The contributions of radial and tangential arcs to the 



distribution of arcs in cluster-centric radius can be clearly 
distinguished. 

• The radial distribution of tangential arcs is very broad 
and extends out to several Einstein radii. Thus, individual 
arcs are not well suited for constraining Einstein radii, in- 
stead, more detailed mass models are needed. 

The work presented here shows that both additional 
structures along the line-of-sight and stellar matter in clus- 
ter galaxies do affect the strong-lensing properties of galaxy 
clusters significantly. Since these effects can boost cluster 
strong-lensing cross-sections by factors > 2, they need to 
be taken into account when using cosmological simulations 
to derive predictions for the giant arc abundance and com- 
paring them to observations. Just by themselves, however, 
they seem to be too small to fully account for the reported 
discrepancy between predicted and obse rved giant arc abun- 
dance in a low-(T8 universe (see, e.g.. iFedeh et al.ll2008l ). 
Thus, other ingredients appear to be needed in addition to 
solve the arc statistics problem. 

An obvious extension to the work presented here is the 
calculation of the expected number of giant arcs in surveys 
with various selection functions. This will require an accu- 
rate knowledge of the source redshift and luminosity distri- 
bution, in particular at high redshifts, as well as a proper 
treatment of magnification bias. 

The discrepancies between our results for the opti- 
cal depths and those based on lower-resolution simulations 
clearly show the need for structure formation simulations 
that accurately resolve the inner regions of clusters on scales 
of a few kpc. On these scales, the complex interplay be- 
tween dark matter, gas, and stars is expected to strongly 
influence the matter distribution. Advances in structure for- 
mation simulations that incorporate physical processes such 
as cooling, star formation, and feedback will help to quantify 
the impact of baryonic physics on the giant arc abundance 
more accurately. 
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